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Abstract 

Above a critical matter density the propagating modes of the neutral kaon system are 
essentially eigenstates of strangeness, but below it they are almost complete eigenstates 
of CP. We estimate the real and imaginary parts of the energies of these modes and their 
mixing at all densities up to nuclear matter density 2 x 10 14 g/cm 3 . In a heavy ion collision 
the strong interactions create eigenstates of strangeness, and these propagate adiabatically 
until the density has fallen to the critical value, whereupon the system undergoes a sudden 
transition to (near) eigenstates of CP. We estimate the critical density to be 20 g/cm 3 , 
and that this density will be reached about 2 x 10 5 fm/c after the end of the collision. 
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Neutral kaon systems are extremely interesting for a variety of reasons [|I[] . When they 
are created in a collision among hadrons or leptons they are essentially in eigenstates of 
strangeness because the strong and electromagnetic interactions conserve flavor. When 
they propagate freely in vacuum the weak interactions are operative and both C and P 
are violated. The eigenstates of the full Hamiltonian are then almost completely flavor 
mixed as short- and long-lived kaons: 



\K S ) = [(l + e)\K°)-(l-e)\K°)\/^2{l + \e\' 

\Kl) = \(l + e)\K°) + (1 - e)\K°) \/J2(l + | e p. (1) 



Here |e| ~ 2 x 10 -3 is the measure of CP violation. When a beam of long-lived kaons is sent 
through ordinary matter, short-lived kaons are generated due to the different interactions 
between the components of the former, namely K° and K°, and atomic nuclei. This 
is called kaon regeneration. The goal of this paper is to study the collective modes of 
propagation of the neutral kaons in dense matter, and to determine the fate of these 
modes after they are created in a collision between large nuclei at high energy. 
We will use a strong interaction basis with 



\K°) = |dS) = |l)=(j) 
\K°) = \ds) = \2)=(° 1 V 



(2) 



Because it is so small and plays no special role in our analysis we shall set e = 0. Then 
\K S ) and \K L ) are eigenstates of CP with eigenvalues +1 and —1, respectively. 

First we do a relativistic analysis. Poles of the propagator determine the time evolu- 
tion of small amplitude excitations. These are obtained from solutions to the equation 

to 2 - k 2 - U vac (uj 2 - k 2 ) - U mat (uj, k) = , (3) 

where the vacuum and matter contributions to the self-energy are indicated. Generally 
it is an excellent approximation to evaluate these 2x2 matrices on the mass shell. This 
means that Il mc is a constant and Tl mat depends on the momentum k only. The usual 
analysis gives 

where A and B are complex numbers. In terms of measurables they are 0: 

A = m K -i(T s + T L )/A 
B = Am/2 + zAr 
m K = ( ms + m L )/2 = 497.67 MeV 
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Am = 3.52 x 1(T 12 MeV 

Ar = r s -r L 

T s = 7.38 x 1(T 12 MeV 
T L = 1.27 x 10" 14 MeV 
Am = (0.478 ± 0.003) Ar. (5) 

The matter contribution is diagonal in flavor and is expressed as 

n mai =(^ °). (6) 

In matter with an excess of baryons over antibaryons the strange and antistrange compo- 
nents, F and F, are different. Of course, this is the origin of kaon regeneration in matter. 
One reason they are different is because the valence antiquark d in the K° can annihilate 
on a valence quark d in the proton or neutron, producing a hyperon. This is not possible 
with a K°. 

Diagonalization of the equation 

( ^ + n)("W(") (7) 



yields the energy eigenvalues 

u± = k 2 + A 2 + B 2 + i(F + F)± ^AA 2 B 2 + |(F - F) 2 (8) 

and eigenstates 

P± = (-X ± V 1 + X 2 ) «± (9) 

where x — (F ~~ F)/AAB. The states evolve in time as exp(— iu±t). In baryon-free 
matter, even at finite temperature, there is no difference between F and F (from the 
point of view of the strong and electromagnetic interactions). Then, with the above 
conventions, f3± = ±a±, and the upper sign corresponds to Kl and the lower sign to K$- 
These are the familiar modes in vacuum except that the square of the mass is shifted by 
F, which generally has both real and imaginary parts. 

At nuclear matter density, 0.155 nucleons/fm 3 , the situation is very different. One 
should expect that \F—F\ is only somewhat smaller than typical nuclear energies, implying 
something on the order of (100 MeV) 2 . Since 

2AB = m K (Am + iAV/2) = (4.19 x 10" 5 MeV) 2 + i(4.28 x 10~ 5 MeV) 2 (10) 

this means that \x\ » 1- In this case, the off-diagonal matrix elements in the self-energy, 
due to the weak interactions, are totally ignorable, and the energy eigenstates are also 
eigenstates of strangeness. 
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Let us now estimate F — F as a function of density. At low to moderate densities 
the self-energy may be expressed in terms of the scattering amplitude for kaons scattering 
on the constituents of the medium, here taken to be an equal mixture of protons and 
neutrons [||]. 



In the first line, np£,{q) is the Fermi-Dirac distribution for a nucleon with momentum q, 
the extra factor of 4 appears because of summation over spin and isospin of the nucleons, 
/ is the forward scattering amplitude in the center-of-momentum frame, and s is the 
usual Mandelstam variable. The angular brackets in the second line denote the particular 
momentum averaging, and is the spatial density of nucleons. If the nucleons are 
not distributed according to a Fermi-Dirac distribution then one ought to use the actual 
momentum distribution; this will affect the numerical value of (/) to some extent. Formula 
(11) is simply a version of the low density virial expansion familiar in statistical physics. Its 
applicability has been discussed in Essentially it requires that the average separation 
of nucleons be larger than both the real part of the forward scattering amplitude and the 
inverse of the average relative momentum between a nucleon and a kaon. One may think 
of this formula quantum mechanically as representing indices of refraction (real part) and 
attentuation (imaginary part). Indeed, it is only a variation of the usual formulae used 
in kaon regeneration sudies. 

There is an identical expression for antikaons with F replacing F and / replacing 
/. The difference / — / has been estimated, or may be inferred, from several sources. 
For example, Eberhard and Uchiyama [|J calculated the real and imaginary parts of the 
forward scattering amplitudes for neutral kaons incident on proton and nuclear targets. 
These calculations were based on measured total cross sections for K + and K~ on protons 
and neutrons, on published values of .ff -nucleon elastic scattering, and charge symmetry. 
From their figures 2 and 3 we estimate that 



Here the averaging is for kaon momenta ranging from 0.3 to 1 GeV, relevant for the 
conditions in hot and dense nuclear matter created in a heavy ion collision. Shuryak 
and Thorsson [||] have calculated the scattering amplitudes for charged kaons incident on 
nucleons based on partial wave analyses of data. From their figures 1 and 2 we estimate 
that 



Here the averaging is for \/s ranging between 1.5 to 1.8 GeV corresponding to the range 
of kaon momenta quoted above. (Unfortunately they have only displayed the real parts 




(11) 



(/-/)« 0.3(1 + z)fm. 



Re</-/)«1.5fm. 
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explicitly.) We do not know the reason for this discrepancy. Fortunately the precise num- 
bers are totally irrelevant for the mixing of the neutral kaons in dense matter, since even 
the smallest estimate of the strong interaction induced self-energy overwhelms that due 
to the weak interactions. For example, using the estimate from Eberhard and Uchiyama, 
we obtain 

F - F « (150 MeV) 2 (l + i) (^-3) , (12) 

which ought to be compared to eq. (10). The fact that the real part of F — F is 
positive implies that the nonrelativistic one-body potential, U = F/2rriK, is relatively 
more repulsive for kaons than for antikaons. This is quite natural for the reasons stated 
earlier. The fact that the imaginary part is positive is also to be expected because the 
imaginary part of the scattering amplitude and cross section are related by the optical 
theorem, a = (Air / k cm )lm f , and antikaons have a much bigger inelastic cross section with 
nucleons than kaons. 

An alternative approach to the kaon self-energy at moderate to high density is to use 
an effective Lagrangian incorporating the relevant degrees of freedom and symmetries |J. 
One finds that 

F(u,k) -- 

F^k) = -^f ps --^. (13) 

J TT J TT 

Here T*kn is the kaon-nucleon sigma term, estimated to be of order 350 MeV. One dis- 
tinquishes the scalar nucleon density, ps = (NN), from the conserved vector density, 
Pn = (Nj°N), although they only begin to differ significantly above twice nuclear den- 
sity. The difference in sign in the above equations results from the fact that the scalar 
density treats particles and antiparticles the same whereas the vector density distinguishes 
particles from antiparticles. Thus 





3p N 


J TT 


4/ 2 ' 




3 pat 


f2 PS 

■1 TT 


4/ 2 



Re (F - F) = = (260 MeV)* [-?-){ ) (14) 




which lies between the estimates from Eberhard-Uchiyama and Shuryak-Thorsson. So 
far, to our knowledge, no one has calculated the imaginary part in this approach. 

What happens to the neutral kaons after production in a heavy ion collision? Based 
on the above analysis of F — F we see that the dimensionless ratio \ is proportional to 
the nucleon density and has a very small imaginary part relative to the positive real part. 
Using the results of Eberhard-Uchiyama we get 

X = 6.3 x 10 12 ( — ^] . (15) 
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From eq. (9) one finds that \a + \ 3> \/3+\ and |/3_| ^> |a_|. Hence the eigenstates are 
|+) = |1) and |— ) = 1 2) with corrections of order 1/x- As the matter expands the density 
decreases until \ becomes very small. Then the eigenstates are the vacuum ones, \Kl) 
and \K S ). //this evolution is adiabatic it would have extremely interesting consequences. 
A K° produced in the collision would evolve into a (rather than into an equal mixture 
of K L and K s ), and a K° would evolve into a K S - Based on valence quark counting U 
one expects that 



Kl 

Ks / observed \ " / in mat ter 





(16) 



That is, the observed ratio of K L to K s should be equal to the ratio of K° to K° produced 
in dense matter, which then is equal to the ratio of charged kaons times the ratio of charged 
pions. The latter product of ratios should not change during the very late dilute expansion 
phase of a heavy ion collision. It has been measured 11 in central Au-Au collisions at the 
AGS at Brookhaven National Laboratory (Eb eam — H GeV/nucleon) with the value 6. 
It has also been measured [[TO) in central Pb-Pb collisions at the SPS at CERN (Ebeam = 
160 GeV/nucleon) with the value 2. Hence, contrary to all other high energy accelerator 
experiments, the ratio of long-lived to short-lived kaons would be far from one! But does 
the neutral kaon system evolve adiabatically? 

To answer this question, consider what happens in the kaon's own frame of reference. 
As time goes on, the surrounding density of matter decreases as 1/t 3 . This naturally 
follows from dimensional analysis, and it also emerges from a calculation of the free 



expansion of an ideal relativistic gas ||TT[. A good estimate of the local nucleon density is 
PN{t) = Po(t /t) 3 . For numerical estimates we shall use p = 0.1 nucleons/fm 3 and t = 10 
fm/c, which are characteristic of central collisions between nuclei of atomic number near 



200 at the AGS and the SPS [12]. Furthermore, the average relative speed in an encounter 
between a kaon and a nucleon will decrease with time. The reason is that the fireball has 
an initial radius of R ~ 10 fm or so. After an elapsed time At = t — to nucleons with 
a relative speed greater than R/ At are unlikely to ever encounter the kaon. If we are 
interested in what happens at late times in the expansion, when the strong and weak 
interactions affecting the neutral kaon system become comparable, the relevant averaged 
scattering amplitudes are not those discussed above. Rather, one can use only the s-wave 
scattering lengths. A compilation of data yields [|13 



fo - f = 0.1 + 0.6i fm , (17) 

with an uncertainty of about 0.1 fm in both the real and imaginary parts. This all leads 
to the difference of one-body potentials being 

U(t) - U(t) = (7.5 + 45z) ( 10 MeV. (18) 
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The problem can be reduced to a two-level Schrodinger equation with a time-dependent 
complex potential. Taking out the kaon mass leads to 

= Hitwt) = ( + Tl) + u ® . ^ Am + ^ Ar - 1 m (19) 

dt n ) flWitJ ^ \Am + \AY -±(T S + T L ) + U(t) ) W " liyj 
The instantaneous energy eigenvalues are 

E ± = \{U + U)- l -(T s + T L ) ± ^(^Am+'-AT^ + 1 -(U-Uy. (20) 

The instantaneous eigenstates are given by eq. (9). The exact solutions can be expanded 
in terms of these, with time-dependent coefficients a±(t), as 

il>(t)=a+(t)\+,t) + a-{t)\-,t) (21) 

where 

i ± ' ( >=(&w) ex p{- i r^> t '}- <22) 

There exists a critical density p c defined by the condition that \x\ = 1. With the above 
input its numerical value is p c = 1.1 x 10 -14 nucleons/fm 3 or about 19 g/cm 3 . For densities 
much greater than p c the strong interactions dominate and the eigenstates of the system 
are eigenstates of strangeness. For densities much less than p c the weak interactions 
dominate and the eigenstates of the system are eigenstates of CP. The critical region may 
be conservatively defined as 8 > \x\ > 1/8 or 8 > Pn/Pc > 1/8. 

The picture that emerges is as follows. For t < 10 5 fm/c the matter expands freely, 
with frequent (on the time scale of 1/Am and 1/AT) interactions of the kaons with the 
nucleons maintaining the kaons in eigenstates of strangeness, namely, K° and K°. For 
10 5 < t < 4 x 10 5 fm/c the system is in a transition region where |%| ~ 1, meaning 
that interactions with the nucleons are comparable in strength with the internal weak 
interactions of the kaons. In this region there is nothing to prevent transitions between 
the states. For t > 4 x 10 5 fm/c the matter is so dilute that the internal weak interactions 
of the kaons dominate and they propagate essentially as in vacuum, namely, as Kl and 
Ks- Thus the states of the neutral kaon system evolve adiabatically at both high and 
low density. The time to pass through the transition region, about 3 x 10 5 fm/c, is so 
short compared to the natural oscillation time of the neutral kaons, 1/Am = 5.6 x 10 13 
fm/c, that this transition may be treated with the sudden approximation. It is at the time 
t ~ 2x 10 5 fm/c and the density p c m 20 g/cm 3 that the strangeness eigenstates K° and K° 
decompose into the eigenstates Kl and Ks of CP. This may be shown mathematically 
from the equations of motion of a±!^\ As a consequence of this sudden transition the 

however, the analysis is more complicated than that given in textbook discussions of time-dependent 
perturbation theory and the adiabatic and sudden approximations because the Hamiltonian is not Her- 
mitian, and the instantaneous eigenstates are not orthogonal in the transition region. 
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observed ratio of long-lived to short-lived kaons should be 1. Finally we should remark 
that the probability for a kaon to decay before the transition region is reached is negligible 
because the lifetime of even Ks is much greater than the time to reach the critical density. 

The phenomenon described here does not happen in elementary particle collisions 
such as e + e~, pp, and pp because the net baryon number is either zero or negligibly small. 
It should be noted that the critical density of 20 g/ cm 3 is characteristic of heavy metalsQ, 
perhaps opening the window on new types of experiments with neutral kaons. 
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